CQ 

C:> 

O 

o 


HDL-CR-  77-020-1 


I 

0 

79 

1 

I 

O 

N> 

O 


Energy  Deposition  Rates  and 
Compton  Electron  Currents  from  Low-Altitude  Bursts 
as  a Function  of  Source  Energy 


I 

I 


cr 

m 

'< 

n 

3 

0 

(8 

X 

3 

n 

0) 

lO 

n 

r“  *< 

0 

o 

Mi* 

s 

o 

o. 

1 

(8 

3>  *D 

CO 

— * 

0 

a 

iA 

t/» 

rt 

n 

C 

X 

Q. 

0 

o 

0 

3 

o 

X 

CD 

;o 

rt 

C 

Ol 

O 

n 

ft 

“1 

i/t 

<8 

rt 

(A 

Qt 

Vi 

D 

01 

a. 

01 

3 

</i 

Cl 

X 

01 

Ol 

o 

n 

rt 

i 

c 

"O 

3 

3 

ft 

O 

O 

rt 

% 

a 

O 

m 

O 

3 

0 

(8 

X 

O 

n 

(8 

rt 

3 

CO 

o 

0 

3 

c 

-1 

r> 

o 

c 

(8 

n 

n 

m 

r8 

3 

3 

(8 

n 

«/) 

(O 

•< 

November  1977 


Prepared  by 


Mathematical  Applications  Group, 
3 Westchester  Plaza 
Elmsford,  N.  Y.  10523 


Under  Contract 

DAAG39-77-C-0020 


Inc. 


D D C 

nr^J 

NOV  3 


1977 


^ UbiL^anns 


O 


U.S.  Army  Materiel  Development 
and  ReodineM  Command 

HARRY  DIAMOND  LABORATORIES 

Adelphi,  Morylond  20783 


FM  ruMK  KlUU^  D«T*MMITnN  UNkMinD 


The  findings  in  this  report  are  not 
to  be  construed  as  an  official  Department 
of  the  Army  position  unless  so  designated 
by  other  authorized  documents. 

Citation  of  manufacturers'  or  trade 
n^unes  does  not  constitute  an  official  in- 
dorsement or  approval  of  the  use  thereof. 

Destroy  this  report  when  it  is  no 
longer  needed.  Do  not  return  it  to  the 
originator . 


Harold  S7/^Schechter y 
Martin  0./Cohen_-^ 


* pcnro^MiHC  oncAMiZATiOM  name  ano  aoomcss 

MATHEMATICAL  APPLICATIONS  GROUP,  INC. 
3 Westchester  Plaza 
Elmsford.  N.  Y.  1052 

cowtaoclTm c omcc  nah#:  ano  aooucss 
Harry  Diamond  Laboratories 

2800  Powder  Mill  Road 
Adel  phi . Md.  20783 


MONtTQAiNC  ACCNCYN  AMC  4 AOOi«E$S<#f  frMi  OtUf) 


1^“"  jmaesfomi:^ 


A17R000^02A1A9 

*i.  TPOaT  OATg  ^ 


6k±>. 


. (mt  thl9  report; 


UNCLASSIFIED 

Tsa  oecl ASSjricArro«''6owNG«AO»N5“ 
SCMCDULC 


n*  oisrnisurioN  statcmcmt  Ui  ihtt 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


(7  OltrniOutiOM  STATCMCmT  fmt  rA«  AAtfrACI  wtiAfAtf  ft*  Bfmck  iP.  ft  BHt9rm*1 


• • tUP^LCMCNTAMY  motes 


HDL  Project  No.  X757E7 


<S  KEY  MOMOS  (Cmt**fftv9  an  rararaA  ft  m*P  fp^tttr  Pf  MacS  waa»»a»> 

Compton  Electron  Currents  SAM-CE  Monte  Carlo  Code 
Energy  Deposition  Radial,  Polar  Components 

Atmosphere  Time  Dependence 

\ Ground  Neutron  Source  Energy  Bands 

AOSTWACt  (C9^**m*m  mm  »aaa»a#  aKa  tt  aacaaaan'  fprtutir  ^ AlacA  nwatAaej 

A series  of  time-dependent  Monte  Carlo  calculations  has  been  per- 
formed to  determine,  as  functions  of  neutron  source  energy  band, 
the  energy  deposition  rates  and  Compton  electron  sources  in  the 
air  due  to  neutron  and  secondary  photon  interactions  from  weapons 
at  1-,  50-,  100-,  and  200-m  heights  of  burst.  The  calculations 
were  made  with  a mooified  version  of  the  SAM-CE  computer  program. 


on 

WW  I JAM  Tl 


lOtTlOM  or  I MOV  It  OMOklTC 


UNCLASSIFIED 


ICCUMTV  CkAMiriCATfON  or  THK  rAOC  tfn  I 


1 


MCu'i’v  CL ic *tion  Of  iMiL  v*ccr>n>««  Om 


^The  energy  deposition  rates  and  Compton  electron  sources  were 
determined  as  functions  of  time  (out  to  100-ms  local  time),  in 
i»9  radial  and  altitude  scoring  bins  surrounding  the  four  iso- 
tropic point  sources.  Answers  were  obtained  for  penetrations 
from  0 to  2. A km  in  the  horizontal  direction  and  from  0 to  1.5  km 
above  the  air-ground  interface. 


The  raw  data  have  been  forwarded  to  HDL  for  smoothing  and  curve 
fitting.  Sample  results  are  presented  and  described  in  this 
report. 


UNCLASSIFIED 


2 


TABLE  OF  CONTENTS 


1.  INTRODUCTION 7 

2.  COMPUTATIONAL  TECHNIQUE  8 


3.  CALCULATIONAL  PROGRAM  11 

3.1  Source  Altitudes  II 

3.2  Source  Energy  Bands ]] 

3.3  Air  and  Ground  Descriptions.  11 

3.4  Cross  Sections  and  Response  Functions 11 

3.5  Special  Scoring  Bins 11 

3.6  Time  Bins 17 

3.7  lmpc"-;nce  Sampling 17 


4.  DESCRIPTION  OF  RESULTS 19 

4.1  Results  Presented 19 

4.2  Variation  of  Time- Integrated  Results 

by  Source  Band  Energy • • 19 

4.3  Time-Dependent  Results  28 


5.  CONCLUSIONS 


62 


Mcasm  Ml 

■ra 

mill  sifiiw^l^ 

m 

latt  SictiH  □ 

■tHOtnca 
M1IFIUII0I .. 

□ 

IT 

iisniiinigi/uiiuiiiin  toia 

tin,  Mi.irtftcuT 


D D C 

NOV  3 1977  I 

JuigSILiJU  IL 

D 


3 


LIST  OF  TABLES 


Table  1 - Source  Energy  Band  Structure,  . 

Table  2 - Composition  Description  ......  I3 

Table  3 ~ Cross  Section  Data  Base ]A 

Table  k - Additional  Subdivision  of  Selected  Regions  Near 

the  Air-Ground  Interface 16 

Table  5 ~ Time  Bin  Structure 18 

Table  6 - Energy  Deposition  in  Region  32  - Weapon  at  200m  HOB 20 

Table  7 ■ Energy  Deposition  in  Region  17  " Weapon  at  200m  HOB 22 

Table  8 - Energy  Deposition  in  Region  53  " Weapon  at  200m  HOB 23 

Table  9 ■ Energy  Deposition  in  Region  32  ■ Weapon  at  Im  HOB 2k 

Table  10-  Energy  Deposition  in  Region  17  ” Weapon  at  Im  HOB 25 

Table  11-  Energy  Deposition  in  Region  53  " Weapon  at  Im  HOB 26 

Table  12-  Gamma  Ray  Energy  Deposition;  Ratio  of 

Im  HOB  to  200m  HOB 27 

k 


LIST  OF  FIGURES 


Figure  1 - Basic  Source-Detector  Geometry  Used 

in  the  Monte  Carlo  Calculations 9 

Figure  2 - Description  of  Scoring  Regions 15 

Figure  3 " Energy  Deposition  Rate  - 50m,  6. 36“^. 07  Mev, 

Neutrons,  Region  I 30 

Figure  4 - Energy  Deposition  Rate  - 50m,  6.36-A.07  Mev, 

Neutrons,  Region  17 31 

Figure  5 “ Energy  Deposition  Rate  - 50m,  6. 36-^. 07  Mev, 

Neutrons,  Region  32 32 

Figure  6 - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

Neutrons,  Region  53 33 

Figure  7 “ Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

Low  Energy  Air,  Region  1 34 

Figure  8 - Enernv  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

Low  Energy  Air,  Region  17 35 

Figure  9 ■ Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

Low  Energy  Air,  Region  32. 36 

Figure  10  - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

Low  Energy  Air,  Region  53 37 

Figure  11  - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

Low  Energy  Ground,  Region  1 38 

Figure  12  - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

Low  Energy  Ground,  Region  I7 39 

Figure  13  “ Energy  Deposition  Rate  - 50m,  6.36-4,07  Mev, 

Low  Energy  Ground,  Region  32 40 

Figure  14  - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

Low  Energy  Ground,  Region  53  ...... .,41 

Figure  15  " Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

High  Energy  Air,  Region  1 42 

Figure  16  - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

High  Energy  Air,  Region  17 43 

Figure  17  " Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev, 

High  Energy  Air,  Region  32 44 


5 


List  of  Tables  (continued) 


Figure  l8  - Energy  Deposition  Rate  - SOm,  6.36-4.07  Mev 
High  Energy  Air,  Region  53  

Figure  19  “ Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev 
High  Energy  Ground,  Region  17 

Figure  20  - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev 
High  Energy  Ground,  Region  1 

Figure  21  - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev 
High  Energy  Ground,  Region  32 

Figure  22  - Energy  Deposition  Rate  - 50m,  6.36-4.07  Mev 
High  Energy  Ground,  Region  53 

Figure  23  - Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Neutrons,  Region  1 

Figure  24  - Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Neutrons,  Region  17 

Figure  25  ” Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Neutrons,  Region  32 

Figure  26  - Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Neutrons,  Region  53 

Figure  27  “ Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Low  Energy  Air,  Region  1 

Figure  28  - Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Low  Energy  Air,  Region  17 

Figure  29  “ Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Low  Energy  Air,  Region  32 

Figure  30  - Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Low  Energy  Air,  Region  53 

Figure  31  " Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Low  Energy  Ground,  Region  1 

Figure  32  “ Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Low  Energy  Ground,  Region  17  • 

Figure  33  " Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Low  Energy  Ground,  Region  32  

Figure  34  - Energy  Deposition  Rate  - 50m,  1.11-0.55  Mev 
Low  Energy  Ground,  Region  53  • 


45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
.60 

.61 


6 


1.  INTRODUCTION 


Mathematical  Applications  Group,  Inc,  (MAGI)  has  performed  a series  of 
time-dependent  Monte  Carlo  calculations  to  determine  the  energy-deposition 
rates  and  Compton  electron  sources  in  air  due  to  neutron  and  secondary  photon 
interactions  from  low-altitude  nuclear  bursts.  These  calculations  were  made 
with  a modified  version  of  the  SAM-CE  computer  program.' 

Isotropic  sources  of  primary  neutron  radiation  were  considered  at  four 
low-altitude  heights  of  burst  (HOB);  1,  50,  100  and  200-m.  Rather  than  treat- 
ing specific  weapon  output  energy  spectra,  as  in  previous  work,^  calculations 
were  performed  for  contiguous  neutron-source  energy  bands  which,  when  combined, 
can  reconstitute  arbitrary  neutron  output  spectra. 


The  energy  deposition  and  Compton  electron  currents  were  determined  as 
functions  of  time  in  49  radial  and  altitude  (R-Z)  scoring  bins  surrounding  a 
burst  point.  _The  atmosphere  was  taken  to  be  homogeneous  at  an  assumed  density 
of  1.11  mg/cm'^.  Answers  were  obtained  for  penetrations  up  to  2.4  km  in  the 
h'  ■'ontal  plane  and  up  to  1.5  km  above  the  ground. 


1.  “.0.  Cohen  aj_. , "SAM-CE:  A Monte  Carlo  Code  for  Three  Dimensional 
Neutron,  Gamma  Ray  and  Electron  Transport  (Revision  5)",  MR-7052-5 
(May  1977). 

2.  M.  0.  Cohen,  H.  S.  Schechter,  and  H.  A.  Steinberg,  "Time-Dependent  Energy 
Deposition  and  Compton  Electron  Currents  from  Three  Selected  Low-Alti- 
tude Bursts",  HOL-CR-76-029-1 /MR-7048  (Aug.  1976). 
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2. 


COMPUTATIONAL  TECHNIQUE 


The  calculations  were  performed  with  a specially  modified  version  of  the 
SAM-CE  Monte  Carlo  program.'  Alterations  to  the  program,  pertinent  to  the 
present  calculations,  include  the  following: 

a.  A special  tracking  procedure  was  used  for  neutrons  with  energies 
low  the  lowest  inelastic  threshold  (for  the  nuclides  which  constitute  the 
and  the  ground).  Also  used  was  a special  thermal  diffusion  model.  These 
procedures  increased  the  computational  efficiency  by  a factor  of  5 to  10. 

b.  Time  dependence  was  recorded  in  local  time  units  defined  as  time  sub- 
sequent to  the  arrival  of  the  uncollided  radiation.  (For  both  primary  neutron 
and  secondary  photon  problems,  local  time  zero  is  defined  as  the  earliest 
possible  arrival  of  photons.) 

c.  Time-dependent  energy  deposition  due  to  neutron  elastic  scattering 
and  photon  Compton  scattering  was  scored  for  all  spatial  regions. 

d.  Compton  electron  sources  were  scored  in  radial  and  polar  bins.  These 
are  now  described: 

Figure  1 shows  the  basic  source-detector  geometry  used  in  the  Monte  Carlo 
calculations. 

A point  isotropic  source  is  located  in  the  Cartesian  coordinate  system  at 
0,0, Z , where  Z is  the  source  altitude,  and  a detector  is  located  at  X^,  Y^, 

Z^.  ^The  X and^Y  axes  define  a plane  parallel  to  the  ground. 

Consider  a vector  score,  F,  at  the  detector  position.  The  score,  in  this 
case  the  average  forward  range  of  a Compton  electron,  can  be  characterized  by 
its  projections  along  the  Cartesian  X,  Y,  and  Z axes.  This  is  not  the  most 
convenient  coordinate  system,  however.  A more_con^enlent  £Oordlnat£  system  is 
defined  by  three  mutually  orthogonal  vectors  I , I , and  lij,  w^ere  1 = the 

radial  vector  which  i£  col  inear  to  the  source-detector  axis;  I - the  poJ_ar 
vector,  where  1^  and  Ip  define  a plane  perpendicular  to  the  ground;  and  = 
the  azimuthal  vector,  which  is  parallel  to  the  ground. 

(It  is  apparent  that  in  a homogeneous  atmosphere,  the  algebraic  sum  of 
all  scores  projected  onto  the  azimuthal  vector  must  vanish.  Hence,  in  the  Monte 
Carlo  calculations,  computer  time  was  not  spent  projecting  individual  scores 
along  this  vector.  The  presence  of  the  ground,  however,  does  produce  net  scores 
along  the  polar  axis,  which  would  otherwise  vanish  in  an  infinite  homogeneous 
a i r med i urn. ) 


be- 
a i r 

ad  hoc 


1.  M.  0.  Cohen  et  al . , "SAM-CE:  A Monte  Carlo  Code  for  Three  Dimensional 
Neutron,  Gamma  Ray  and  Electron  Transport  (Revision  5)",  MR-7052-5 
(May  1977). 
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At  each  photon  Compton  scattering  event  in  the  Monte  Carlo  game,  special- 
ly prepo.-'d  built-in  data  tables  were  used  to  generate  the  average  forward 
range  of  the  resulting  recoil  electron.  (These  tables  were  obtained  by  folding 
in  the  known  recoil  electron  angle-energy  probability  distribution  as  a function 
of  electron  energy  and  then  by  resolving  the  results  along  an  axis  parallel 
to  the  incident  photon.) 


Subsequently,  the  average  forward  range  was  resolved  along  the  radial  and 
polar  scoring  axes  described  above. 


e.  The  secondary  gamma  problems  were  divided  into  four  component  calcu- 
lations. Photon  energy  deposition  and  Compton  sources  were  tallied  separately 
for : 


(1)  Photons  generated  by  the  interactions  of  "high  energy" 
neutrons  (E>0.1  MeV)  in  the  air. 

(2)  Photons  generated  by  the  interaction  of  h^gh  energy 
neutrons  in  the  ground. 

Photons  generated  by  the  interaction  of  "low  '•nergy" 
neutrons  (E<0,1  Mev)  in  the  air. 

(4)  Photons  generated  by  the  interaction  of  low  energy 
neutrons  in  the  ground. 


The  reason  for  this  separation  of  secondary  photons  was  to  simplify  sub- 
sequent smoothing  and  curve-fitting  analyses  to  be  performed  by  HDL. 

f.  In  order  to  further  simplify  the  analysis  task  by  HDL,  all  scores 
were  written  on  magnetic  tape  (in  a fixed  format)  and  delivered  directly  to  HDL. 

g.  For  all  combinations  of  HOB  and  source  energy  band,  computer-generated 
plots  of  secondary  photon  energy  deposition  rates  were  provided  for  selected  re- 
presentative spatial  regions. 
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3.  CALCULATIONAL  PROGRAM 


I 


3. 1 Source  Al t i tudes 

Separate  series  of  calculations  were  performed  for  HOB  values  of  1,  50, 

100  and  200  meters."  The  results  obtained  (see  below)  indicate  that  inter- 
polation procedures  will  be  adequate  to  scale  the  results  to  any  other  HOB  in 
the  ground-to-200  meter  interval. 

3.2  Source  Energy  Bands 

The  source  energy  band  structure  is  given  in  Table  1.  The  100m  HOB  prob- 
lem was  run  first,  and  subsequent  analyses  of  the  sensitivity  of  the  results 
to  source  energy  showed  that  the  final  band  structure  (as  used  for  the  1,  50, 
and  200m  HOB  runs)  would  be  adequate. 

3.3  Air  and  Ground  Descriptions 

The  air  composition,  as  used  in  the  calculations,  is  given  in  Table  2.  It 
represents  dry  air  at  a density  of  1.11  mg/cm^,  corresponding  to  a height  above 
sea  le-'cl  of  approximately  900  m.  A homogeneous  model  of  the  atmosphere  was 
assumed . 

The  ground  composition  which  was  used  is  also  given  in  Table  2.  It  corres- 
ponds to  dry  Nevada  test  soil  at  a density  of  1.7  g/cm^. 

No  attempts  were  made  to  assess  the  effects  upon  the  results  of  moist  air 
or  ground  compositions  other  than  dry  Nevada  test  soil. 

3.^  Cross  Sections  and  Response  Functions 

The  cross  section  data  base  for  the  computations  is  given  in  Table  3.  The 
average  range  data  for  Compton  electrons,  referred  to  in  Section  2,  were  devel- 
oped by  the  authors  after  consultation  with  various  knowledgeable  experts  in 
the  field. 


3.5  Special  Scoring  Bins 

The  atmosphere  was  divided  into  stacked  and  concentric  cylinders  for  the 
purposes  of  scoring  and  importance  sampling  (see  Section  3.7).  The  spatial 
mesh  was  300  m (33.3  g/cm^).  A finer  subdivision  of  the  atmosphere,  however, 
was  used  closer  to  the  source. 

Figure  2 shows  the  spacing  of  the  scoring  regions  and  the  identifying  re- 
gion number  associated  with  each  scoring  bln.  (The  unlabeled  regions  were  pre- 
sent in  the  monte  Carlo  calculations  for  backscatter Ing  purposes  only.) 

'n  order  to  investigate  the  behavior  of  the  scored  quantities  near  the  air- 
ground  interface,  regions  2,  i*,  10,  and  12  were  further  subdivided  as  shown  in 
Tabic  If. 


* 


The  1 m HOB  is  used  to  represent  a ground  burst. 
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TABLE  1 


Source  Energy  Band  Structure 


**Brackets  show  multiple  bands  which  were  subsequently  collapsed  to 
a single  band. 
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TABLE  2 


Composition  Descriptions 


AIR  (dry) 


Photon  Transport 

N;  3.2501X10"^  atom/g.  air  XIO^^ 
0:  8.7201X10"^  '' 

Ail:  1.9l«36X10"^  " 


Neutron  Transport 
N:  3.2655X10"^ 

0:  8.7610X10"^ 


GROUND 

(dry  Nevada  test 

soil) 

Photon 

and  Neutron  Transport 

H: 

5.7'*8X10‘^ 

atom/g.  ground  XI 0 

0: 

2.046X10"^ 

M 

Si : 

6.822X10”^ 

II 

AH  : 

2.872X10"^ 

II 

Dens i ty 

Air: 

Ground: 

1.11  mg/ctj? 

1 .70  g/cm^ 
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TABLE  3 


Cross  Section  Data  Base 

•I 


Neutron  Transport  and 

Photon 

Element 

Photon  Production  Data 

Transport  Data 

Hydrogen 

DNA  l4ll*8-Mod  2 

ENDF/B  as 
disseminated  by 

RSIC  as  Tape 

N i t rogen 

DNA  1*1 33-Mod  k 

DLC-7  D 

Oxygen 

DNA1»  134-Mod  2 

A1 uminum 

DNA  4135-Mod  2 

(all  elements 
including  Argon) 

Sill  con 

DNA  4l5l-Mod  2 

Source  on  vertical  Axis 
at  1 , 50,  100  m (Reg.  1 ) 
or  200  m (Reg.  5 } 


FIGURE  2 - Description  of  Scoring  Regions 
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NOMINAL  REGION 

SUBDIVISIONS 

ALTITUDE  RANGE 

(ref«r  to  Fig.  2) 

(Region  Number) 

(Meters) 

2 

so 

0-  25 

51 

25-  50 

52 

50-100 

2 

100-150 

k 

53 

0-  25 

5^ 

25-  50 

55 

50-100 

100-150 

10 

56 

O 

1 

vn 

57 

25-  50 

58 

50-100 

59 

100-150 

10 

150-300 

12 

60 

0-  25 

61 

25-  50 

62 

50-100 

63 

100-150 

12 

150-300 
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3.6  Time  Bins 

The  primary  neutron  and  secondary  photon  scores  were  obtained  in  local 
time  bins  out  to  a 100  millisecond  time  cutoff.  The  local  time  structure  is 
given  in  Table  5. 

3 . 7 Importc ' -e  Sampl i ng 

The  probler s , be  solved  involved  difficult  Monte  Carlo  importance  samp- 
ling situations.  .nswers  were  required,  with  low  statistical  uncertainties, 

over  many  regions  of  a large  mul t i -d imens ioned  phase  space.  Some  of  the  tech- 
niques used  to  ensure  adequate  solutions  included  the  following: 

(a)  Special  tracking  for  neutrons  below  the  inelastic  threshold  and  a 
special  thermal  diffusion  model  - as  described  in  Section  2,  item  (a). 

(b)  Energy  importance  sampling  to  discriminate  against  low  energy  neut- 
ron collisions  which  would  be  unlikely  to  generate  secondary  photons. 

(c)  Spatial  importance  sampling  to  generate  a sufficient  number  of  neutron 
high  energy  interactions  at  shallow  penetrations  (since  the  photons  which  they 
generaLw  dominate  the  early  temporal  ranges). 

(d)  Spatial  '-portance  sampling  to  obtain  adequate  solutions  at  the 
deep  penetrations  and  to  "push"  a sufficient  number  ot  neutron  events  towards 
the  ground  and  towards  the  special  subdivided  volumes  of  Table  A. 

(e)  Directional  importance  sampling  to  discriminate  against  ground- 
generated photons  initially  headed  downwards  deeper  into  the  ground. 
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TABLE  5 


Time  Bin  Structure 


Time  Bin 

Time  Interval 

1 

( 0 

- 0. 

2 

( 0.100 

- 0. 

3 

( 0.215 

- 0. 

4 

( 0.464 

- 1. 

5 

( 1.00 

- 2. 

6 

{ 2.15 

- 4. 

7 

( 4.64 

- 10. 

8 

(10.0 

- 21. 

9 

(21.5 

- 46. 

10 

(46.4 

-100. 

1 1 

( 0.100 

- 0. 

12 

( 0.215 

- 0. 

13 

( 0.464 

- 1. 

14 

( 1.00 

- 2. 

15 

( 2.15 

- 4. 

16 

( 4.64 

- 10. 

17 

(10.0 

- 21. 

18 

(21.5 

- 46. 

19 

(46.4 

-100. 

*Read:  0.100  x 10  ^ to  0.215  x 10  ^ sec 
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(sec) 
100)  E“6 
215)  E-6* 
46A)  E-6 
00  ) E-6 
15  ) E-6 
64  ) E-6 
0 ) E-6 

5 ) E-6 

4 ) E-6 

0 ) E-6 

215)  E-3 
464)  E-3 
00  ) E-3 
15  ) e-3 
64  ) E-3 
0 ) E-3 

5 ) E-3 
4 ) E-3 
0 ) E-3 


DESCRIPTION  OF  RESULTS 


A. 


^ . 1 "esults  Presented 

The  final  results  presented  in  this  report  will  be  sharply  limited  for 
the  following  reasons: 

(a)  The  primary  objective  of  the  studies  was  to  determine  energy  deposi 
tion  rates  and  Compton  current  sources  in  the  proscribed  phase  space  bins  and 
to  forward  these  raw  data  to  HDL  for  subsequent  smoothing  and  curve-fitting 
analyses.  This  objective  was  met  and  results  for  each  source  band,  at  each 
HOB,  have  been  delivered  by  MAGI  to  HDL  as  formatted  output  on  magnetic  tape. 

(b)  Much  of  the  qualitative  analyses  for  the  previous  study,  involving 
fission  and  lA  MeV  source,'  apply  to  the  band-source  results  of  this  study 


as  well. 

These  include: 

(1) 

general  temporal  shapes  of  the  neutron  and 
components  of  the  energy  deposition  curves 

four  photon 

(2) 

general  temporal  shapes  of  the  four  photon  components 
of  the  radial  and  polar  Compton  current  sources 

(3) 

air-ground  interface  effects 

(M 

range  and  HOB  effects  on  each  of  the  above. 

Hence,  for  such  analyses,  the  reader  is  referred  to  the  previous  work.  However, 
in  the  sections  which  follow,  some  of  the  salient  new  results  are  presented. 

k.2  Variation  of  Time-Integrated  Results  by  Source  Band  Energy 

Table  6 displays  the  time-integrated  energy  deposition  from  the  neutron  and 
the  four  photon  components  in  Region  32  (see  Figure  2)  as  a function  of  source 
energy  band.  The  results  are  for  the  200-m  HOB,  integrated  in  time  out  to  the 
100  ms  cutoff.  Region  32  was  selected  for  this  first  set  of  presented  results 
because  it  is  at  an  intermediate  distance  from  the  source,  and  is  600  to  900-m 
above  the  air-ground  interface. 

Table  6 shows  that  the  high  energy  - air  photon  component  is  the  most  im- 
portant contribution  to  the  total  dose  for  the  higher  energy  source  bands.  How- 
ever, as  the  neutron  source  energy  drops  below  the  air  inelastic  thresholds, 
this  component  completely  disappears. 


M.  0.  Cohen,  H.  S.  Schechter,  and  H.  A.  Steinberg,  "Time- Dependent  Energy 
Deposition  and  Compton  Electron  Currents  from  Three  Selected  Low-Altitude 
Bursts",  HDL-CR-76-029-1/MR-701*8  (Aug.  1976). 
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ENERGY  DEPOSITION  IN  REGION  32 
Weapon  at  200-m  HOB 

(Integrated  from  local  time  0 to  100  mil 


The  high  energy  - ground  photon  component  is  less  important  than  the 
high  energy  - air  component  at  the  highest  energies  but  extends  down  to  lower 
neutron  source  energies  due  to  the  lower  inelastic  thresholds  in  the  ground. 

The  two  low-energy  photon  components  are  smaller  than  the  high  energy  - 
air  photon  component  at  the  highest  neutron  source  energies.  However,  they 
remain  markedly  constant  in  their  contribution  (per  source  neutron)  all  the  way 
down  to  the  lowest  neutron  source  band.  In  general,  the  total  contribution 
from  the  low  energy  - air  component  is  more  than  twice  that  of  the  low  energy  - 
ground.  This  would  not  have  been  our  conclusion  if  the  time  cutoff  had  not  been 
extended  out  to  100  ms,  in  contrast  with  the  earlier  work  which  was  out  to  10  ms 


Table  6 also  shows  that  neutron  energy  deposition  is  more  than  S0%  of  the 
total  deposition  in  the  neutron-source  energy  range  of  2.35  to  6.36  MeV.  Above 
6.36  MeV,  neutrons  produce  many  high-energy  air  gamma  rays,  by  inelastic  scatter 
ing,  which  dominate  the  deposition  of  energy.  Below  2.35  MeV  the  neutrons  do 
not  have  much  energy  themselves  but  do  produce  capture  gamma  rays  which  domin- 
ate the  deposition  of  energy.  (The  trends  noted  in  this  paragraph  hold  for  all 
heights  of  burst  and  all  spatial  regions.  For  example,  see  Tables  7 through  11 
which  fol low. ) 

Table  7 shows  conditions  similar  to  those  of  Table  6,  except  that  the  spa- 
tial region  is  17,  which  is  closer  to  the  source  (and  to  the  ground)  than  is 
region  32.  Similar  trends  are  seen.  It  is  noticed,  however,  that  the  low 
energy  - ground  component  has  been  enhanced,  vis-a-vis  the  low  energy  - air 
component  and  now  exceeds  the  latter  for  some  of  the  source  energy  bands. 

Also  note  that  since  the  region  is  closer  to  the  source,  neutron-energy  deposi- 
tion dominates  over  a wider  range  of  source  energy  bands. 

Table  8 shows  conditions  for  region  53  which  is  one  of  the  subdivisions  of 
region  k (see  Table  k)  and  which  extends  down  to  the  ground.  Trends  similar  to 
those  noted  above  are  again  observable.  Due  to  their  nearness  to  the  interface, 
both  photon  ground  components  have  been  enhanced  vis-a-vis  the  photon  air  com- 
ponents. Indeed,  the  low  energy  - ground  component  clearly  dominates  the  low 
energy  - air  component. 


For  example,  at  the  100-m  HOB  an  analysis  of  the  low  energy  - air  component 
for  energy  band  no.  ^ (''-fission  source)  and  band  no.  13  (''-1^  Hev  source) 
show  65%  and  59^,  respectively,  of  the  t Ime- Integrated  dose  in  the  extra  10 
to  100  ms  time.  By  contrast,  the  values  are  26%  and  18%,  respectively,  for 
the  low  energy  air  - ground  component. 


ENERGY  DEPOSITION  IN  REGION  I7 
Weapon  at  200-m  HOB 


Read:  O.8I  x 10  ev/ (cm^* source  neutron) 


ENERGY  DEPOSITION  IN  REGION  53 
Weapon  at  200um  HOB 

(Integrated  from  local  time  0 to  100  milliseconds) 


;ead:  0.26  x 10  ev/(cm  'source  neutron 


ENERGY  DEPOSITION  IN  REGION  53 


Read:  0,23  x 10  ev/ source  neutron) 


TABLE  12 

GAMMA  RAY  ENERGY  DEPOSITION;  RATIO  OF  1-m  HOB  to  200-m  HOB 


COMPONENT 


ENERGY 

BAND 
(Mev  ) 

High 

Energy-Ai r 
Region 

17  32  53 

High 

Energy-Ground 

Region 

17  32  53  1 

Low 

Energy-Air 

Region 

17  32  53 

Low 

Energy-Ground 

Region 

17  32  53 

.0335-  .11 

g 

D 

3.0 

2.4 

n 

.11  - .55 

2.2 

2.2 

.55  - 1.11 

0.8 

0.5 

B 

1.9 

1.9 

0.3 

1.11  - 1.83 

0.8 

0.5 

1.1 

2.4 

1.7  1 

0.7 

1.83  - 2.35 

8.1 

7.2 

3.'» 

1.0 

0.8 

1.0 

2.0 

1.3 

0.5 

2.35  - '♦.07 

0.5 

0.1 

1.2 

13. 

8.0 

1.0 

0.8 

0.4 

2.3 

1.3 

1.3 

^4.07  - 6.36 

O.l^ 

0.6 

0.8 

6.3 

11.3 

1.2 

1.3 

0.7 

2.3 

1.8 

0.6 

6.36  - 8.19 

0.5 

0.5 

0.6 

9.5 

8.5 

0.8 

0.8 

0.7 

2.1 

1.3 

1 .0 

8.19  -15.0 

O.A 

0.5 

0.7 

9.0 

6.'* 

0.8 

lA 

0.5 

2.7 

2.0 

1.1 

Tables  9,  10  and  11  are  repeats  of  Tables  6,  7 and  8,  except  that  the 
height  of  burst  has  been  lowered  from  200-m  to  1-m.  For  regions  17  and  32, 
the  ground-phcton  component  is  seen  to  be  greatly  enhanced.  (Also  see  Table  12, 
below).  For  region  53,  the  1 i ne-of-s 1 ght  from  the  source  to  the  detector  for 
the  1-m  case  is  so  close  to  the  horizontal  that  the  ground  component  may  actually 
be  decreased  in  some  cases  (also  see  Table  12,  which  followS)  vis-a-vis  the  200-m 
HOB  case. 


Table  12  displays  the  ratio  of  the  gamma  ray  deposition  components  for 
the  1-m  and  200-m  heights  of  burst.  For  regions  17  and  32,  it  is  seen  that  lower- 
ing the  source  height  from  200-m  to  1-m: 

(a)  reduces  the  high  energy  - air  component  by  about  one-half, 

(b)  greatly  enhances  the  high  energy  - ground  component  by 
about  a factor  of  ten, 

(c)  reduces  somewhat  (sometimes  less  than  by  one-half)  the  low 
energy  - air  component,  and 

(d)  enhances  the  low  energy  - ground  component  by  about  a factor 
of  2 which  is  a smaller  enhancement  than  for  the  high 
energy  - ground  component. 


Region  53  is  an  exception,  as  explained  in  the  preceding  paragraph.  The 
results  for  a 1-m  HOB  and  a 200-m  HOB  are  much  closer  to  each  other  than  was  the 
case  for  regions  17  and  32. 

i».3  Time-Dependent  Results 

Figures  3 through  which  follow,  display  representative  energy  deposition 
rates  in  regions  1,  17,  32,  and  53  ^or  the  50-m  height  of  burst.  (Region  1 con- 
tains the  source.)  Trends  similar  to  those  of  the  50-m  HOB  were  found  for  the 
other  three  altitudes. 

In  all  figures,  the  ordinate  values  are  in  units  of  eV/ (cm^* sec* source 
neutron),  and  the  abscissa  values  are  in  units  of  seconds  (of  local  time)'.  Ex- 
cept for  the  neutron  deposition  histograms,  the  symbols  on  the  plots  represent 
the  calculat ional  statistical  accuracy  to  the  nearest  tens  of  percent.  For  ex- 
ample, a "3"  represents  30^  statistics.  For  neutrons,  where  only  zeros  (0)  appear, 
statistics  are  not  available. 


Figures  3 through  22  display  the  results  for  the  A. 07  through  6.36  MeV  range 
which  is  representative  of  the  higher  neutron  source  energies. 


Figures  3 through  6 show  the  neutron  results,  for  this  source  energy  band, 
for  regions  1,  17,  32,  and  53,  respectively.  Figure  3,  for  region  1,  shows  the 
energy  deposition  existing  as  local  time  goes  to^zero.  Jhis  is  to  be  expected, 

-f • « 1 


since  region  1 contains  the  source.  Between  10 
a rapid  decay  with  time  and  disappears,  for  all 


and  10  5 seconds  the  pul^se  begins 
practical  purposes,  by  10  ^ 


Recall  that  local  time  zero  is  defined  as  the  earliest  possible  time  of  arrival 
of  secondary  gamma  radiation. 
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seconds.  Figures  4-6  show  the  arrival  of  the  neutron  pulse  at  'vlO  seconds  and 
then  similar  decay  schemes.  In  fact,  the  pulse  shapes  of  Figure  3“6  are  quali- 
tatively similar  to  the  corresponding  curves  for  other  heights-of-burst,  spatial 
regions,  and  source  bands. 

Figures  7"10  show  the  low  energy-air  results  for  regions  1,-  17,  32  and  53. 
respectively.  Although  slight  variations  exist,  these  curves  are  representative 
of  the  results  found  for  all  heights-of-burst,  spatial  regions  and  source  bands. 
The  curves  rise  steeply  at  about  10  3 seconds,  may  show  a slight  rise  for  a short 
time  period,  (e.g.,  see  Figure  8)  and  then  decay  away  rather  slowly  in  time  out 
to  the  cutoff  of  this  problem  which  is  10  ' seconds. 


Figures  11-14  show  the  results  for  the  same  four  spatial  regions  for  the 
low-energy  ground  component.  Neutrons  slow  down  much  more  rapidly  in  the  ground 
than  in  the  air  and  thus  they  are  a I so^captured  more  quickly.  Hence,  the  low- 
energy  ground  component  arrives  at  10  to  10  5 seconds,  the  exact  value  de- 
pending upon  the  HOB  and_the  altitude  of  the  detector  region.  Also  note  that 
by  the  cutoff  time  of  10  ' seconds,  the  magnitude  of  this  component  has  been 
reduced  from  the  peak  value  by  several  orders  of  magnitude.  The  decay  of  the 
low  energy  - ground  component,  as  shown  in  Figures  11-14,  is  characteristic  of 
all  heights-of-burst,  spatial  regions,  and  source  energies. 


Figures  15~io  display  the  high  energy-air  component  for  the  four  spatial 
regions.  Since  neutrons  can  collide  at  time  zero  (0)  and  produce  high  energy 
gamma  rays  which  fly  uncollided  out  to  a scoring  region,  this  component  of  the 
energy  deposition  is  seen  to  exist  at  the  earliest  local  times.  After  remaining 
relatively  constant  for  several  orders  of  magnitude  of  local  time  (for  which 
theoretical  arguments,  not  given  here^i-can  be_provided)  this  component  decays 
away  rapidly  and  vanishes  at  about  10  ^ to  10  **  seconds.  The  shape  of  these 
figures  is  characteristic  for  all  heights  of  burst  and  all  spatial  regions.  How- 
ever, they  are  very  much  a function  of  source  band  energy.  Indeed,  the  compon- 
ent does  not  even  exist  below  the  inelastic  threshold  in  air  (e.g.,  see  Table  6). 


Figures  19“22  display  the  high  energy  - ground  component  for  the  four 
spatial  regions.  The  shape  of  the  pulse  is  characteristic  of  all  heights  of 
burst  and  spatial  regions,  but,  as  was  the  case  for  the  high  energy  - ground  com- 
ponent, a function  of  the  source  and  energy.  Note  that  the  pulse  is  of  rela- 
tively short  duration,  rising  at  about  KT®  seconds  (thus  being  a function  of 
the  HOB,  i.e.,  the  flight  distance  from  source  to  the  ground)  and  decaying  away 
by  about  10  **  seconds. 


Figures  23-34  are  for  0.55“I.ll  MeV , which  is  one  of  the  low-energy  source 
bands.  They  correspond  to  Figures  3“14,  respectively,  for  the  4.07-6.36  Mev 
band.  There  are  no  corresponding  figures  to  Figures  15-22  since,  as  noted  above, 
the  high  energy  photon  components  do  not  exist  for  this  source  energy  band. 

Comparisons  of  Figures  23-34  with  the  corresponding  Figures  3-14  bear  out 
the  comments  made  above  that  the  characteristic  shapes  of  the  neutron  and  low 
energy  photon  components  are  relatively  independent  of  neutron  source  energy 
band. 
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5.  CONCLUSIONS 


The  ilr-ry  Diamond  Laboratories  (HDL)  required,  by  neutron  source  energy 
band,  the  energy  deposition  rates  and  Compton  currents  due  to  bursts  at  I,  50, 
100  and  200  «>.  The  spatial  range  was  out  to  2. A km  in  the  horizontal  direction 
and  from  0 to  1.5  km  above  the  air-ground  interface.  The  temporal  range  was 
out  to  100  ms.  Through  the  careful  use  of  importance  sampling  and  some  special 
coding,  MAGI  has  used  the  SAM-CE  Monte  Carlo  program  to  generate  the  raw  data, 
with  sufficient  accuracy,  in  the  important  phase  space  regions.  These  data 
have  been  forwarded  to  HDL  for  subsequent  processing  (see  below).  Typical  raw 
data  results  have  been  presented  (for  energy  deposition)  in  Section  k to  high- 
light the  more  important  trends. 

The  HDL  in-house  program  smooths  and  establishes  analytical  fits  to  the 
raw  data  supplied  by  MAGI.  To  this  end,  the  four-way  split  in  the  secondary 
gamma-ray  sources  not  only  helped  to  understand  the  nature  of  the  transport 
processes  which  are  involved,  but  will  also  make  the  smoothing  and  curve  fitting 
considerably  easier  to  perform.  The  fact  that  the  characteristic  shapes  of 
the  temporal  curves  are  relatively  independent  of  height  of  burst  and  of  spatial 
region  (and,  in  some  cases,  of  the  source  energy)  will  also  make  the  task  of 
HDL  notably  easier  to  perform. 
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